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Abstract. The use of ultra-thin, i.e., monolayer films plays an important role for the 
emerging field of nano-fluidics. Since the dynamics of such films is governed by the 
interplay between substrate-fluid and fluid-fluid interactions, the transport of matter 
in nanoscale devices may be eventually efficiently controlled by substrate engineering. 
For such films, the dynamics is expected to be captured by two-dimensional lattice-gas 
models with interacting particles. Using a lattice gas model and the non-linear diffusion 
equation derived from the microscopic dynamics in the continuum limit, we study two 
problems of relevance in the context of nano-fluidics. The first one is the case in which 
along the spreading direction of a monolayer a mesoscopic-sized obstacle is present, 
with a particular focus on the relaxation of the fluid density profile upon encountering 
and passing the obstacle. The second one is the mixing of two monolayers of different 
particle species which spread side by side following the merger of two chemical lanes, 
here defined as domains of high affinity for fluid adsorption surrounded by domains of 
low affinity for fluid adsorption. 



PACS numbers: 68.15.+e, 68.43.Jk 
1. Introduction 

There are substantial efforts to miniaturize chemical processes by using microfiuidic 
systems. The "lab on a chip concept" integrates a great variety of chemical and physical 
processes into a single device in a similar way as an integrated circuit incorporates 
many electronic devices into a single chip [Tj. These microfiuidic devices do not only 
allow for cheap mass production but they can operate with much smaller quantities 
of reactants and reaction products than standard laboratory equipments. This is 
particularly important for solutions containing rare and expensive substances, such as 
certain biological materials, and for toxic or explosive components [Sj. Even though most 
microfiuidic devices available today have micron sized channels, further miniaturization 
is leading towards the nano-scale IS]. Besides meeting technical challenges, new 
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theoretical concepts are needed to understand the basic physical processes underlying 
this new technology jU Ej- Whereas the ultimate limits for miniaturization of 
electronic devices are set by quantum fluctuations, in a chemical chip these limits 
are determined by thermal fluctuations and can be explored by methods of classical 
statistical mechanics. 

At the sub-micron thickness scale, recent experiments of liquid spreading on 
atomically smooth surfaces [7J |H1 E] , performed with volumes of the order of nano-liters, 
have clearly shown by means of dynamic ellipsometry or X-ray reflectivity measurements 
that precursor films with molecular thickness and macroscopic extent advance in front 
of the macroscopic liquid wedge of the spreading drop. ( Thin, i.e., of the order of 100 
nm, precursor films spreading ahead of the macroscopic droplet have also been observed 
experimentally |10j.) The occurrence of molecularly thin precursor films with a similar 
spreading dynamics has been evidenced very recently also for immiscible metal systems 
in three regimes: solid drops with a solid film, solid drops with a liquid film, and liquid 
drops with a liquid film [TTJ [TJj. Theoretical work (see Refs. [El CU 113 HE] and 
references therein) combined with an impressive number of Molecular Dynamics (MD) 
and Monte Carlo (MC) simulations (see Ref. ^1 and references therein) addressed 
the mechanisms behind the extraction and the experimentally observed t 1 ^ 2 asymptotic 
time dependence of the linear extent of the precursor films on chemically homogeneous 
substrates H31 E3 EE] ■ This led to a good understanding of the spreading dynamics 
and of the intrinsic morphology of the films. Based on these results, more complicated 
issues can be addressed such as, e.g., the spreading behavior of monolayers exposed to 
chemically patterned substrates ^3 1201 , or the question of mixing of different fluids at 
the nano-scale that we shall present below. 

The organization of the paper is as follows. In Sec. El we briefly present the 
lattice gas model of interacting particles and discuss the rules defining the microscopic 
dynamics and the nonlinear diffusion equation derived from it in the continuum limit. 
Section |3] is devoted to a qualitative discussion of the results obtained for the case of 
monolayer spreading in the presence of mesoscopic obstacles, with a particular focus 
on the relaxation of the density profile upon encountering and passing the obstacle. In 
Sec. 0] we discuss the mixing of two species during spreading of monolayers following 
the merger of two chemical lanes, and we conclude with a brief summary of the results 
in Sec. 

2. Fluid monolayers on homogeneous substrates 

Recently, we have studied the structure of a monolayer, which extracts from a reservoir 
[TH] and spreads on a flat, chemically homogeneous substrate, by using a lattice gas 
model of interacting particles as proposed in Refs. [131 HHj. Since we shall use this 
model as a starting point for our present study, for clarity and further reference we briefly 
describe the defining rules of the model and the non-linear diffusion equation obtained 
from the microscopic dynamics within the continuum limit. A thorough analysis of this 
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model is presented in Ref. [T5] . 

(a) We choose a homogeneous substrate such that the spreading occurs in the x — y 
plane. The half-plane x < is occupied by a reservoir of particles at fixed chemical 
potential which maintains at its contact line with the substrate — positioned at the line 
x = — an average density Co (defined as the number of particles per unit length in 
the transversal y direction). At time t — 0, the half-plane x > is empty. There is no 
imposed flow of particles from the reservoir pushing the extracting film. 

(b) The substrate-fluid interaction is modeled as a periodic potential forming a lattice 
of potential wells with coordination number z (z = 4 for a square lattice) and lattice 
constant a. The particle motion proceeds via activated jumps between nearest-neighbor 
wells; evaporation from the substrate is not allowed. The activation barrier Ua 
determines the jumping rate Cl = vq exp[— C/a/^bT], where v is an attempt frequency 
defining the time unit, ks is the Boltzmann constant, and T is the temperature. 

(c) The pair interaction between fluid particles at distance r is taken to be hard-core 
repulsive at short range, preventing double occupancy of the wells, and attractive at 
long range, —Uo/r 6 for r > 1, resembling a Lennard- Jones type interaction potential. 
Here and in the following all distances are measured in units of the lattice constant a 
and therefore are dimensionless. The selection of the nearest-neighbor well into which 
a particle attempts to jump, i.e., the probability p(r — > r'; t) that a jump from location 
r will be directed toward the location r', is biased by the fluid-fluid energy landscape 
and is given by 

exp{l[U(r;t)-U(r';t)}} 
P(r-r;*) = ^ , (1) 

where Z(r;t) = exp|^[?7(r; t) — U(r'; £)] j is the normalization constant and 

r' ,\r'— r|=l 

l/P = k B T, 

0(r;t) = -U (2) 

r',0<|r'-r|<3 

and rj(r';t) G {0,1} is the occupation number of the well at r' at the time t. The 
summation in Eq. (J2J has been restricted to three lattice units for computational 
convenience. This corresponds to the cut-off generally used in Molecular Dynamics 
simulations for algebraically decaying Lennard- Jones pair-potentials. The rates 

^r^r';t = &p(r — > r' \ t) (3) 

for the transitions from r to neighboring sites r' satisfy 

uj r ^ rl -t = CI. (4) 

r',\r'— r|=l 

Thus for any given particle at any location the total rate of leaving a potential well is 
determined only by the fluid-solid interaction characterized by Ua, is time- independent, 
and equals Cl. 
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Neglecting all spatial and temporal correlations, i.e., assuming that averages of 
products of occupation numbers 77(7*; t) are equal to the corresponding products of 
averaged occupation numbers p(r;t) = (i](r;t)), where (. . .) denotes the average with 
respect to the corresponding probability distribution P({p(r; t)}) of a configuration 
{t](r;t)}, one can formulate a mean-field master equation for the local occupational 
probability, i.e., the number density p(r;t) |18j . In the continuum limit of space and 
time (At — > 0, a —>■ 0, Q^ 1 — > 0, D = Qa 2 /4 finite) for the master equation, by 
taking Taylor expansions for p(r — > r') and p(r'; t) around r and keeping terms up to 
second-order spatial derivatives of the density p(r;t) |H1|2T|, one obtains the following 
nonlinear and nonlocal equation for p(r\t) fZ*2\ 12 3j : 

d tP = D V [Vp + /3p(l- p)VU] + 0(a 2 ) (5) 



where 



U(r;t) = (U(r;t)} = -U ^ (6) 

r",0<|r"-r|<3 ' ' 

is replacing t/(r; t) in the definition ((TJ for p(r — > r'). 

Being nonlinear and, due to the term involving the interaction potential U(r;t), 
nonlocal, Eq. (|SJ) cannot be solved analytically and in most of the cases even the 
computation of a numerical solution is a difficult task. However, assuming that the 
density p(r; t) is a slowly varying function of the spatial coordinates the potential U(r; t) 
may be expanded as 

r',0<\r'-r\<3 

^-U p(r;t) T^h^ + ° {a2) > (?) 

r',0<|r'-r|<3 ' ' 

which leads to the local equation 

d tP = A)V{[1 - g W p(l - p)} Vp} + 0(a 2 ), (8) 
where Wq = (3Uq, and g = l r l~ 6 is a geometrical factor depending on the lattice 

l<\r\<r c 

type (e.g., square, triangular, etc.) and on the cut-off range of the potential. For the 
present case of a square lattice and a cut-off at r c = 3 one has g ~ 4.64. 

Rescaling time as t — > r = Dot and defining an effective diffusion coefficient 

D e (p) = l-gW p(l-p), (9) 

Eq. (jHJ) may be written in the usual form of a diffusion equation: 

d T p = V [D e (p)Vp] + 0(a 2 ) . (10) 

The functional form of D e (p) (Eq. Q) implies that for Wo > 4/g there will be values 
Pi of the density for which D e (pi) < 0. For parameters such that Wq < 4/g, Eq. (fTDj) 
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is a proper diffusion equation (though non-linear), while for Wq > 4/g instabilities are 
expected in the range of densities for which D e (pi) < 0, i.e., for p { 6 (p~ , p+) where 



It is known |2"2"1 |2U 123 12S| that these instabilities lead to discontinuities in the density 
profile ("shocks"), i.e., they correspond to the formation of sharp interfaces. For the 
model defined by the rules (a)-(d), the value for the threshold interaction strength for 
which such interfaces emerge is predicted by the continuum theory as W — 4/g ~ 0.86, 
which is significantly smaller than the lower bound estimate W > 1 from KMC 
simulations. We attribute this to the mean-field character of the derivation of the 
continuum equation. Therefore it is necessary to include particle-particle correlations 
into the mean-field description. Since the dynamics is possible only by jumps into empty 
sites, one can argue that for z = 4 the summation in g should include at most three 
contributions from nearest neighbor sites. This leads to g ~ 3.64 and an estimate for 
the threshold interaction ~ 1.1, in good agreement with the KMC results. For the 
rest of the analysis we shall use this corrected value of g. Additional support for this 
corrected value is provided by the analysis of the density profiles [T5] . 

The constraint, that the reservoir keeps the mean density at x = at a fixed value 
Co, implies the boundary condition 

p(x = 0,y;t) = C . (12) 

Depending on the particular system under study, additional boundary conditions may 
have to be satisfied. For example, for the case studied in Ref. [IE], in the absence of 
formation of interfaces, i.e., for interactions Wq < W^p and for large times, the density 
on the advancing edge X(t) could be considered as fixed and equal to Ci, 

p(x = X(t),y;t) = C 1 , (13) 

where C\ = 0.11 as inferred from the kinetic Monte Carlo (KMC) simulations. (This 
boundary condition (Eq. p3)l) naturally occurred also in the theory of Burlatsky et al 
[To].) For that system, the absence of boundaries along the y— direction and the y- 
independence of the boundary conditions at x = and x = X(t) leads to an effectively 
one-dimensional problem and to a scaling solution p(x, t) — C(X = xj y/i). The analysis 
of Eq. (Jini) depends on whether W < or W > W • As shown in Ref. jig], 

in both cases the solutions are in excellent agreement with those obtained from KMC 
simulations; typical results are shown in Fig. ^ While for the liquid-on-solid systems 
mentioned in the Introduction these intrinsic density profiles have not been measured 
yet, data of such density profiles are available for the immiscible metal systems studied 
in Ref. E] and they are in at least good qualitative agreement with the theoretical 
ones. Since the present model appears to provide a simple but realistic description of 
a fluid monolayer spreading on a homogeneous substrate, it is natural to use it as a 
starting point to address more complex problems, such as the spreading of monolayers 
on designed chemically heterogeneous substrates, or the mixing of monolayers. 
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Figure 1. Asymptotic scaling solution p(x,t) = C(A) for Wo = 1.0 (regular solution) 
and (see the inset) Wq = 1.4 (shock solution) with A = x/\fDot. The theoretical results 
are obtained from Eq. I|10f) (solid lines). The KMC results (open circles) correspond 
to the time t — 2 x 10 6 (in units of Vq ) which is close to the asymptotic limit 
(see Ref. |18j). A s denotes the position of the discontinuity as obtained by analytical 
theory [IB] - 



3. Diffusive spreading around mesoscopic obstacles 

In order to apply the model described in Sec. [21 to study the spreading of a monolayer 
around a mesoscopic-size obstacle, we add to the rules (a)-(d) in Sec. II the following 
ones: 

(e) The obstacle is taken to be a square-shaped domain T> of side length h centered at 
(x — d > h/2, y = 0) [see also Fig. [2(a)]. This domain is composed of sites with very low 
affinity for the fluid particles. The activation barrier Uv for jumps from sites outside T> 
to those inside T> is taken to be much larger than Ua, such that the boundary dT> of T> 
acts effectively hard wall. 

(f) A sink reservoir occupies the region x > L, where L ^> 1 and L ^> d + h, and 
maintains at its contact line with the substrate, positioned at the line x = L, an average 
density (number of particles per unit length in the transversal y direction) C\ = 0. 

Under these assumptions, the density profile p(x, y, t) as the solution of Eq. (|T0|) 
fulfills the initial condition 

p(x,y,0) = C o e(-x), (14) 

where Q(x) denotes the Heaviside step function, and the boundary conditions 

p(0,y,t) =C , 

p{L,y,t) =0, (15) 
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Figure 2. (a) Schematic drawing of a substrate with high affinity for fluid adsorption 
patterned with a square domain T> [length h, centered at (x — d > h/2, y = 0)] 
of very low affinity for fluid adsorption obstructing the spreading. The substrate is 
in contact with a reservoir of particles located at x < from which a monolayer is 
extracted and spreads along the x direction. A sink of particles is located at x = L. 
There are periodic boundary conditions in the y— direction. (b),(c) Density profiles 
in the vicinity of the obstacle from numerical integration of Eq. llUf) with initial and 
boundary conditions given by Eqs. (fHI) - (fR))) . respectively, for a spreading monolayer 
whose edge just encounters an obstacle (r = 200) (b) and has just passed the obstacle 
(t = 2000)(c). The parameters used in the simulations are Co = 1.0, Wo = 0.9, d = 20, 
h = 10, L x = 200, L y — 30. The right box shows the color coding for the density. 



jn \&D = 0. 

The current j is given by [see Eq. (jlOj) ] 

j = -D e (p)Vp. (16) 

Note that in a numerical study the system necessarily has a finite size L y along 
the y direction. We will use periodic boundary conditions and sizes L y ^> 1 such that 
h/Ly < 1 but not negligible (mesoscopic-size obstacle) and L y — h^> r c , such that the 
boundary dT> of the obstacle is sufficiently far away from the edge of the simulation box 
to avoid finite-size effects. 

In Figs. Efb) and (c) we present typical results for the density profiles in the vicinity 
of the obstacle obtained from the numerical integration of Eq. (jlOJ) with initial and 
boundary conditions given by Eqs. (|nj) - (|T6|) . respectively, for a spreading monolayer 
whose edge just encounters an obstacle (r = 200) (b) and has just passed the obstacle 
(r = 2000) (c). Several conclusions can be drawn from visually inspecting Figs. E{b) and 
(c). Upon approaching the obstacle, the boundary condition of zero normal current at 
the boundary of the obstacle (reflecting wall) leads to an increase in the density in a 
region at the front of the spreading monolayer, as shown by the forward bending of the 
iso-density lines [see, e.g., the yellow band in Fig. I2^b)]. Upon passing the front-edge 
of the obstacle, the iso-density lines become straight, and once they reach the end of 
the obstructed region they bend again, this time backwards. This indicates that upon 
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passing the obstacle the spreading tends to proceed faster in the regions far from the 
obstacle, while at the obstacle the iso-density lines are pinned until they cover the whole 
edge on the back of the obstacle [see, e.g., the light blue band neighboring the green 
region in Fig. I2K C )]- Once this is realized, the iso-density line detaches from the back- 
edge of the obstacle, and the bending slowly relaxes [see, e.g., the boundary between the 
light and dark blue regions in the top region of Fig. 121(c)], while the spreading continues; 
far away from the obstacle, the iso-density lines become again straight. 

We end this section by noting that the maximum linear extent of the region where 
the iso-density lines are deformed, as well as the survival time of these deformations can 
be used as quantitative measures to describe the relaxation of the density perturbations 
induced by the obstacle as a function of the inter-particle attractive interactions, as well 
as of the scaled size h/L y of the mesoscopic obstacle (assuming that this is the most 
relevant geometrical parameter). The results of this analysis will be the subject of a 
forthcoming paper. 

4. Diffusive mixing of two fluid monolayers composed of different species 

The model described in Sec. El can be used to study the mixing of two spreading fluid 
monolayers composed of different species A and B, if one assumes that the two species 
interact with the substrate in such a way that the same lattice structure of potential 
wells, eventually with different depths (i.e., different escape rates f2j), can accommodate 
both types of species. Assuming a square lattice of lattice constant a and assuming 
the on-site hard core repulsion between any two particles such that double occupancy 
remains forbidden, the equations satisfied by the densities pj(x,t), where j 6 {A, B}, 
are obtained from Eq. (jHJ) by replacing U with Uu + Uij, where Uu is the potential due to 
same-species interactions while Uy, j ^ i, is the potential due to interactions between 
different species, Dq with D{ = f2ja 2 /4, and changing the single-occupancy term from 
1 - p(r, t) to 1 - p A (r, t) - p B (r, t): 

d tPi = AV [V Pl + P Pi (1 - Pi - Pj )V(U u + U i:j ] ,i,je {A, B}, j^i. (17) 

We note here that there is no summation over the same indices, and we also note that 
in the case of identical species, i.e., Da = D b and U aa = Ubb = Uab, by adding 
Eq. (jl7j) for i = A and the corresponding one for i = B one finds that as expected the 
density p = Pa + Pb satisfies Eq. (0). Assuming that the long-ranged parts of the pair- 
interactions A — A, B — B, and A — B are of the same Lennard- Jones form (see Sec. El) 
only with different strengths Uq , i, j G {A, B} } one can repeat the same argument as 
that following Eq. (JHJ) to reduce Eq. fTTj) to a local one, 

d tPi = AV [Vp, - g Pi (1 - Pi - Pj ) {WuVpi + WijVpj] , i, j E {A, B}, j^i, (18) 

where the notation Wij = (3Uq^ has been introduced. 

We consider a T-junction patterned onto a planar, rectangular substrate of size 
L x x 2L y as described below [see also Fig.|2Ka)]. The spreading monolayers are extracted 
from reservoirs of particles A and B, respectively, which maintain constant line densities 
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Figure 3. (a) Schematic drawing of a T-shaped pattern with high affinity for a fluid 
in contact with two reservoirs for particles of species A and B, located at y = —L v 
and y = L y , respectively. The two monolayers spread initially in opposite directions 
along the y axis, and upon meeting they start to mix and to invade the empty stripe 
of width h located at x > d, \y\ < h/2. The rectangular domains T>\ and V 2 consist 
of sites with very low affinity for fluid adsorption. A sink of particles is located at 
x = L x . (b), (c) Results for the mixing parameter ca{t,t) near the T-junction 
from numerical integration of Eq. i|18[l with initial and boundary conditions given 
by Eqs. I|20 [l -I|21 [) for the cases of attractive inter-species interaction Wab = 0.6 (b) 
and repulsive inter-species interaction Wab = —0.6 (c), respectively, at r = fO 3 , for 
Waa = Wbb = 0.7, and Qf = = 1.0. (d), (e) Results for the density p A (r,r) of 
A particles for attractive inter-species interaction Wab = 0.6 (d) and repulsive inter- 
species interaction Wab = —0.6 (e). The other parameters are the same as in (b) and 
(c). Note that there is no color-coding for Pa{t, t). 



Cq ,b at the lines C\ = (y = —L y , < x < d) and C 2 = (y = L y , < x < d), respectively. 
The domains £>i = {(x,y)\y < —h/2 A x > d} and V 2 = {(x,y)\y > h/2 A x > d} 
represent sites with very low affinity for the fluid particles of either type, such that 
similarly to the situation in Sec. El the boundaries of these domains effectively act as 
hard walls, confining the spreading onto the two lanes forming the inverted T-junction. 
Finally, we assume that at the foot x = L x , \y\ < h/2 of the T-junction there is a sink 
for particles of both species. Under these assumptions, the functions pA,B(x,y,t) as 
solutions of Eq. (j!8j) fulfill the initial condition 

p A (x,y,0) = Cq, (x,y) G Ci , 

p B (x,y,0) = C*,(x,y)eC 2 , (19) 
Pa,b(x, y, 0) = 0, otherwise, 
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and the boundary conditions 

PA(x,y,t) \ Cl = C, A 



o • 

B 



p B (x,y,t)\ C2 =C^ (20) 
p A (L x ,y,t) = p B (L X} y } t) = 0, 



3n' B \dVi,2 — 



'1,2 
:A I -B 



where the current j (j is obtained by exchanging the labels A <-> B) is now given by 
[see Eq. (JTHD] 

j A = Vp A - 5 p A (1 - p A - p B ){W AA Vp A + V^ B Vp B ) . (21) 

In the following we focus on the effect of the A-B interaction on the dynamics of 
mixing of otherwise identical monolayers, i.e., we choose W AA = Wbb, D a = D B , and 
Cq = Cq ; the results discussed in the following correspond to the particular choice 
Cq = Cq = 1, while the parameter D A is absorbed into the variable r = D A t. The 
geometrical parameters are fixed to L x = 500, L y = 20, and d = h = 20. The mixing 
will be characterized by the ratio 

ex(r,r)= . >f' T \ . (22) 

(with the convention = if p A + p^ = 0), which is close to 1 in A- rich regions, close 
to zero in B-rich regions, and close to 1/2 in regions where mixing is accomplished (i.e., 
Pa ^ Pb ^ 0). 

In Figs. Efb) and (c) we present typical results for c A (r,r) from numerical 
integration of the coupled set of equations given by Eq. (fTHJ) with the initial and 
boundary conditions Eqs. (j2U|l - (j2T|l for the case of attractive interactions W AA = Wbb = 
0.7 and attractive inter-species interaction W AB = 0.6 (b), respectively repulsive inter- 
species interaction W AB = —0.6 (c). From these figures it is clear that within the stripe 
forming the leg of the T-junction there is almost perfect mixing. 

Surprisingly, at first glance the result seems to be almost independent of the sign 
of the A-B interaction, and there is only a weak dependence on the strength W AB 
of this interaction, except for the extension of the A-rich and B-rich regions near the 
corners of the T-junction. However, this behavior can be easily rationalized in view 
of the fact that, as shown in Ref. [IB], the structure and dynamics of the spreading 
of a one-component monolayer of particles with inter-particle attraction W AA = 0.7 
is well described also by the "effective boundary force" theory of Burlatsky et al [15J, 
which disregards interactions within the bulk of the monolayer. The reason for this is 
that the density in the spreading monolayer is relatively low, except for the region near 
the reservoir, and thus the system is too dilute to be influenced by the inter-particle 
attraction. Therefore, in the initial stages of spreading and mixing, the stripe is invaded 
by low density phases of A and B which mix independently of their mutual interaction. 
This scenario is well supported by the comparative analysis of the corresponding density 
profiles p A (r,r) in the cases W AB = 0.6 and W AB = —0.6, respectively, shown in 
Figs. Efd) and (e): almost everywhere in the stripe the density of A particles is low, 
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in the range of to 0.25, and thus the mixture behaves as a dilute, non-interacting 
two-dimensional gas. Finally, we note a weak dependence of the extension of the A-rich 
and B-rich regions near the corners of the T-junction on the sign of the Wab interaction 
[compare Figs. El (b) and (c), respectively (d) and (e)]. 

The above scenario holds for all values Waa = Wbb < 0.9 and |Wab| < Waa 
that we have tested, and thus the (tentative) conclusion is that the symmetric T- 
geometry would ensure practically perfect mixing (but without any possibility to control 
the spatial extent or the spatial distribution of mixing) for two-component monolayers 
with a similarity of the interactions between like species. However, further calculations 
should be carried out before definite conclusions can be drawn concering this issue. For 
example, for the case Waa, Wbb > 11 we expect sharp interfaces to emerge in each 
of the two spreading monolayers, which might eliminate the "mixing through the low 
density front" mechanism discussed above. Simulations for these ranges of values for 
the inter-particle interactions turned out to be extremely time consuming, and work is 
still in progress to elucidate this point. 

5. Summary 

A lattice gas model of interacting particles and the corresponding nonlinear diffusion 
equation derived from its microscopic dynamics in the continuum limit provide a simple 
but realistic description of fluid monolayer spreading on a homogeneous substrate. Based 
on previous results for spreading on a homogeneous substrate, here we have extended 
this model to address two more complex problems: the spreading of monolayers around 
obstacles (Sec. EJ) and the mixing of monolayers (Sec. HJ). These are simple examples of 
spreading on chemically designed substrates. 

For the case of monolayer spreading in the presence of a mesoscopic obstacle, 
the results obtained from the numerical integration of the nonlinear diffusion equation 
(Eq. (fTUJ) ) with initial and boundary conditions given by Eqs. (jHJ)-([TBJ) show that the iso- 
density lines are bent and pinned by the obstacle during the spreading of the monolayer 
around it. For a fixed geometry, a fixed density of the reservoir, and fixed substrate-fluid 
and inter-particle interactions the spatial and temporal extent of this bending in front 
of and behind the obstacle can be used as measures for the relaxation of the density 
profile upon passing around the obstacle. 

As an example for the mixing of two species in the course of spreading of two 
monolayers at the merger of two chemical lanes, we have discussed the case of a T- 
junction geometry. We have focused on the effect of the A-B inter-species interaction on 
the dynamics of mixing of otherwise identical monolayers. Surprisingly, so far our results 
lead to the conclusion that the symmetric T-geometry together with the similarity of the 
same-species interaction ensures practically perfect mixing (but without any possibility 
to control the spatial extent or the spatial distribution of mixing) for two monolayers of 
different species, independently of the sign or the strength of the A-B interaction. Only 
the extension of the A-rich and B-rich regions near the corners of the T-junction exhibits 
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differences. This behavior reflects the fact that in the initial stages of spreading and 
mixing the stripe is invaded by low density phases of A and B, which mix independently 
of the inter-species interaction; the resulting mixture is a dilute, quasi non-interacting 
two-dimensional gas. Because this complete mixing has occurred at early stages, when 
repulsion does not play a role, and since the continuum equation does not contain any 
noise terms, demixing or segregation is not observed in the present calculations, although 
it is expected to occur for strongly repulsive A-B interactions. 
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